Many body exchange effects close to the s-wave Feshbach resonance in two-component 

Fermi systems: is a triplet superfluid possible? 
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We suggest that the exchange fluctuations close to a Feshbach resonance in a two-component Fermi 
gas can result in an effective p-wave attractive interaction. On the BCS side of a Feshbach resonance, 
the magnitude of this effective interaction is comparable to the s-wave interaction, therefore leading 
to a possible spin-triplet superfluid in the range of temperatures of actual experiments. We also 
show that the particle-hole exchange fluctuations introduce an effective scattering length which does 
not diverge, as the standard mean-field one does. Finally, using the effective interaction quantities 
we are able to model the molecular binding energy on the BEC side of the resonance. 
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In the atomic Fermi gases such as 40 K and 6 Li, the use 
of Feshbach resonances has opened the possibility of ex- 
ploring the very interesting limit for which the mean-field 
approximation predicts a smooth crossover from BEC 
to BCS pairing as one goes through the resonance. At 
low energies, the inter-atomic interaction is very well de- 
scribed by the s-wave scattering length, a s . Moreover, no 
direct interactions are possible in the triplet channel. In 
fact, higher-order expansions in the scattering length are 
suppressed at very low temperatures and the symmetry 
of the wave function, due to Pauli exclusion, does not 
allow s-wave scattering for fermionic atoms in the same 
spin channel [l|]. 

Although the scattering length in the two-body prob- 
lem is diverging, it is instructive to consider the possi- 
bility of pairing in the higher-order scattering channels 
due to exchange fluctuations. It is also not clear whether 
atomic systems behave as Fermi liquids (FL), or how sim- 
ilar they are with high T c superconductors (HTSC) or 
any other strongly correlated systems. 

In this Letter, we want to show two things. Firstly, 
that it is possible to build a Fermi liquid theory (FLT) 
in the atomic Fermi gases, particularly in the BCS re- 
gion. This formalism explains the basic features of these 
gases like the scattering lengths, and possibly, the bind- 
ing energies in the strongly interacting regime, which is 
not accessible by simple perturbation theory. Secondly, 
we show that important contributions can arise in higher- 
order momentum channels. The resulting triplet pairing 
is comparable, in magnitude, to the s-wave one, and the 
correspondent triplet superfluid transition temperature 
is within experimental reach. Since the triplet interac- 
tion can only occur through the fluctuations induced by 
the strong interactions with the other spin channel, we 



focus our discussion primarily on the induced term. Its 
contribution may give an instability of the Fermi sea for 
(quasi-)particles with equal spins, and leads to a possible 
transition to a triplet superfluid. We will keep the discus- 
sion quite general, since our approach may be of interest 
to other fields. Then, we specify to the cold atom physics 
case, as we proceed. 

In a Fermi liquid at sufficiently low temperatures (T c < 
T < Tp, where T c is the BCS transition temperature 
and Tp is the Fermi temperature), it was shown in Q 
that it is possible to separate the Fermi liquid parame- 
ters, /pp', which by construction do not contain any zero 
sound terms Q , into two sets of terms in the limit when 
p rsj p' (here, we assume the general notation q = (q, lo) 
and q = |q|). One term is the direct interaction of the 
quasiparticles (QP), and the other is the crossed term of 
the particle-hole contribution. Dropping for the moment 
the spin indices, the Fermi liquid parameters are 
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where Z Pi are the residues of the single particle Green's 
functions at the pole of the QP, ?? p is the Fermi distri- 
bution function, U is some interaction, and Sn pi /5U P2 
is related to the response function and can be obtained 
from the QP transport equation 0]. Restoring the spin 
indices, we denote the first term by d^, and the second 
term by I™, . The many-body effects in the QP interac- 
tion are therefore separated into two contributions: 
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where d 1 ^,, the direct term, includes only the diagrams 
which are not particle-hole reducible and is equivalent 
to the T-matrix in the particle-particle channel. The in- 
duced term, I™, , has contributions from the exchange 
of virtual collective excitations among the quasiparticles, 
i.e. density, spin-density, current, spin-current fluctua- 
tions to name a few (see [2J for full details) . The implicit 
assumption, as it is for all Fermi liquid theories, is that 
all the relevant processes occur on the Fermi surface. 

Consider, now, any quantity, say /|' a . This is related 
to its counterpart f™, by the definition 
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where F^' a = 7V(0)/|' a , N(0) is the density of states at 
the Fermi surface. The superscript s (a) indicates the 
symmetric (anti-symmetric) contribution with respect to 
the spin, the subscript £ indicates the Legendre compo- 
nent, and Pi is the Legendre polynomial . Then, by 
expanding the Bethe-Salpeter equation for the QP inter- 
actions in the limit w/|q| — > in a rotationally invariant 
system into Legendre polynomials, it can be shown that 
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Here, A s ^ = N(0)a'^' i> is the symmetric (anti- 
symmetric) Legendre components of the scattering am- 
plitudes of the quasiparticles. Note that these scattering 
amplitudes differ from the bare scattering amplitudes, 
since they contain the many-body effects of the theory 
through the QP interactions /. Given that Eq. ^ is a 
non-perturbative result, it remains valid even when Ft 
diverges, since the Ag remain finite. The only approxi- 
mation at this point has been in assuming a Fermi liquid 
and the low energy and momenta limits. From Q, it 
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FIG. 1: Schematic diagrammatic relation between the Lan- 
dau parameters and the scattering amplitudes defined in 
Eqs. (|2I4I5I6|1 . The first term on the RHS of / is the direct in- 
teraction d, and the second is the induced term /. Notice that 
when neglecting the induced terms, this description reduces 
to the RPA. 



follows that 
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where q' 2 = |p— p'| 2 = fc|(l — cos6*l) and cos#l = PF>' is 
the Landau angle, and Xo(q') is the density-density cor- 
relation (Lindhard) functions (see Q). Including I > 1 
terms is straightforward in this model, but only leads 
to small corrections to the results. For the direct term 
D in the low temperature limit, the particle-particle 
T-matrix is proportional to the bare s-wave scattering 
length a s . Since D is then angle independent, it con- 
tributes only in the £ = momentum channel, given 



by D s 



-N(0)U/2, where U = 4irha s /m is 



some on-site interaction. The direct interaction is anti- 
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symmetric and obeys the Pauli principle D = 0. We 
purposely neglect the remaining diagrams in the particle- 
particle channel, since we are mostly interested in the 
induced interaction driven by the exchange of collective 
excitations between the quasi-particles. In fact, at these 
temperatures, this is the only way a triplet interaction 
can arise in the same spin channel. We observe that 
the form of the direct term depends on the model used, 
whereas the induced term does not. Still, we emphasize 
that the resulting scattering lengths, calculated through 
the Bethe-Salpeter equation, have the correct symme- 
tries and conserve the Pauli principle through the Lan- 
dau sum rule, given by ^2g{A^ + Af) = 0, which is not the 
case for the random phase approximation (RPA). Look- 
ing at the diagrams in Fig. [T] may help understand the 
differences. The RPA lack of the exchange terms in the 
particle-hole channel implies an inconsistent treatment of 
the QP interactions, since the scattering amplitudes are 
not properly anti-symmetrized. The consequences of this 
will appear clear below. 

We now apply the above Fermi-liquid formalism to the 
specific case of the cold atomic Fermi gases, in particular 
to 40 K. In these systems, the scattering length can be 
varied by tuning the system close to a magnetic Feshbach 
resonance The s-wave scattering length a s (denoted 
in the figures as ag are ) varies as 
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where Bo denotes the magnetic field value of the Fesh- 
bach resonance, AB is the width of the resonance, and 
dbg is the background scattering length. Since most 
of the experimental systems deal with broad resonances 
only, the contribution of the molecules from the closed 
channel can be neglected 0]. 

We solve self-consistently Eqs. (|5I6|) by varying the di- 
rect interaction a s and use Eq. ^ to obtain the scat- 
tering amplitudes A)' 3 " for I — 0,1. We then use these 
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FIG. 2: S-wave scattering lengths (units of Bohr radii) as a 
function of the magnetic field on both sides of the Feshbach 
resonance Bo, in 40 K, using data from The particular 

scattering length is denoted in the legend. The 50% error bars 
in the experimental data have been removed for clarity. We 
also plot the RPA scattering length for comparison. 



scattering amplitudes to construct the singlet and triplet 
airing amplitudes in the well-known s-p approximation 
and call this the effective singlet scattering length 
a° ff . In Fig. [2] we show the results for the bare (Eq. [7], 
dashed line) and effective (thick line) s-wave scattering 
lengths calculated in this model for 40 K on both sides of 
the Feshbach resonance B = Bq. The most important 
feature is that the effective scattering length does not 
diverge as the resonance is approached, while a s does 
diverge. Our results are very close in slope and magni- 
tude to the experimental values. On the other end, far 
from Bo, the effective and mean- field scattering lengths 
are comparable. We note that the presence of a strong 
p-wave interaction would influence the background scat- 
tering length abg- However, since in that channel there is 
no resonance, the many-body effects will give a negligible 
contribution and we can safely assume the background 
scattering length to be constant. 

Also shown in Fig. [5] is the result predicted by RPA 
(thin line), which clearly fails to capture the correct 
physics as the resonance is approached. The divergence 
of the RPA scattering length implies the emergence of 
two new ground states on cither side of the Feshbach 
resonance. On the BEC side, this would correspond to 
the Stoner instability or the onset of ferromagnetism. On 
the BCS side, this would correspond to phase separation. 
We note that neither of these two instabilities has been 
observed experimentally. 

We now turn to the BCS side of the resonance where 
the scattering lengths are negative, and compare the 
effective triplet (a° s ) and effective singlet scattering 
lengths scaled by the Fermi momentum kp . For magnetic 
fields near the resonance, the strength of the triplet and 
the singlet potentials are actually comparable in mag- 
nitude, as shown in Fig. [3J Therefore, although direct 
(triplet) pairing in the same hyperfine state is suppressed 



initially at these low temperatures, the exchange of col- 
lective excitations upon approaching the resonance drives 
a substantial attractive interaction in the triplet channel 
through the induced interactions. Note that there is al- 
ways attraction in the triplet channel on both sides of 
the resonance, independent of the sign of the bare inter- 
action U. At low enough temperatures, it is not obvious 
that one can disregard the possibility of having a triplet 
superfluid near resonance (on the BCS side). In fact, be- 
low the triplet transition point, it seems quite reasonable 
that these two many body states will compete, as Fig. [3] 
suggests. 

It is important to compute the critical temperatures 
expected for the various pairing instabilities. The expres- 
sion is very similar to the BCS one (see 0] and Q). The 
singlet transition temperatures with our effective scatter- 
ing amplitude could be as large as T c smg ~ .7Tp, while 
the triplet T* rip ~ .2 Tp, if we use Tp as the cut-off scale. 
Note that these critical temperatures are quite high and 
that this is due to the use of the high-energy cut off. 
Also, there are numerous indications that singlet tran- 
sition temperatures arc of the order of .2Tp. This in- 
troduces a proportionality relation which gives a (upper 
bound of the) triplet transition temperature, within our 
framework, of ~ .05Tp. These temperatures are already 
obtainable in current experiments. We also mention that 
in the limit when the partial waves get very large, i.e., 
far away from the resonance, our approach gives the 
Gorkov and Melik-Barkhudarov critical temperature 0, 
since the particle-hole corrections become unimportant 
in this regime. 

Lastly, we should remark that the present calculations 
hold for both equal populations or for very small polar- 
izations, m = {n) — n^)/n, where n^^' is the majority 
(minority) particle density, and n is the total density. For 
m <C n, the corrections to the Fermi liquid parameters, 
which are quadratic in m, are, in fact, negligible. Re- 
cent experiments [h]] have opened up the possibility of 




FIG. 3: Scattering amplitudes on the (attractive) BCS side of 
the resonance Bo. In our model, the effective (unitless) s-wave 
singlet fcFof = A^ s and triplet kFal s = A% s amplitudes are 
finite and of similar magnitude at the resonance, while the 
bare interaction fefls = A^ arc diverges, according to Eq. Q. 
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FIG. 4: Binding energies on the BEC side in 40 K, com- 
pared with data taken from[TlI[. The density used here is 
n = 5.8 x 10 13 cm -3 , with AB = 9.7 G and the Feshbach 
resonance, marked above with the vertical line, occurs at 
Bo = 224.21 G. Ef is the binding energy calculated from 
the effective scattering length in our theory, and Eh is the 
bare (mean-field) binding energy. 
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FIG. 5: Gap squared using the weakly interacting limit equa- 
tions as derived in but replacing the effective mass and 
effective interaction from our theory to the bare ones. U is 
the bare interaction strength as defined in the text. V is re- 
lated to the gap squared as in [(|, although our units do not 
correspond to those of the reference. 



exploring the triplet interactions of the system. Since the 
singlet BCS state is still stronger than the triplet one, in 
order to see the triplet transition, it is probably necessary 
to suppress the singlet superfluid. Indeed, even a small 
polarizations, at low enough temperatures, might cre- 
ate the possibility of a triplet superfluid in that channel. 
More likely, the presence of an external field can establish 
a preferential direction and favor the triplet pairing, sim- 
ilar to the Ai phase in 3 He. We emphasize, though, that 
the triplet pairing is not only interesting in the superfluid 
phase, but also in the normal one, as it will contribute to 
the properties of the system. In fact, its thermodynamic 
properties, which we will discuss somewhere else, can be 
largely affected and therefore provide the experimental 
tools to verify the Fermi liquid behavior of the system 
close to resonance. 

Up to this point, we have assumed that the correc- 
tions due to the particle-particle contribution are not 
relevant. In the normal phase, this assumption is plau- 
sible on the BCS sideQ, but on the BEC side can be 
justified only very close to the resonance and/or for tem- 
peratures k^T > E\, (the experimental data are taken 
at T/Tp ~ .4, although in this paper we are consider- 
ing only the corrections due to the quantum fluctuations 
T < Tp ), where E\> is the binding energy. Thus, deep 
into the BEC regime, our theory breaks down. On the 
other hand, we recover the bare scattering length as soon 
as we get far away from the resonance and therefore any 
FLT assumption is irrelevant. 

Thus, we might expect the theory to give a rough esti- 
mate of the scattering lengths in the intermediate region 
as well. We therefore compute the binding energy of the 
bound state in the open channel. It is calculated using 
the standard mean-field formula E^ = —h 2 / ma 2 , where 
a s is the bare s-wave scattering length. Since we lack 



a better estimate of the corrections to this formula due 
to the many-body effects, we simply replace a s with a° s 
and m with m* = 1 + -Ff/3, the effective mass in FL 
theory. The results arc shown in Fig. 21 The agreement 
with the experimental data is quite surprising, but can 
be explained in terms of an effective Hamiltonian, which, 
in the spirit of Landau's theory, progressively transforms 
the bare particles into quasiparticles and the bare scatter- 
ing length becomes the effective one as the interaction in- 
creases. In this sense, one can adopt the same mean-field 
formula of the binding energy, and indeed, the effective 
binding energy profile reduces to the mean-field one in the 
weakly interacting regime. We also note that the Fermi 
liquid parameters represent only a part of the mean-field 
shift, since they do not contain the zero sound channel 
contribution. The full contribution is instead given by 
the full effective scattering length. 



Engclbrccht et aZ.[l2| calculated the energy gap equa- 
tion in the weakly interacting limit. They correctly 
pointed out that a comparison of the gap with the full 
solution should show roughly the extent of the strongly 
interacting regime. It is also interesting to note that the 
BEC-BCS crossover behavior is lost by using the weak- 
limit gap solution. In Fig. (|5|) we plot the energy gap cor- 
responding to the weakly interacting BEC and BCS lim- 
its, but with our effective scattering lengths and masses. 
It is clearly seen that the crossover is re-established in 
terms of the quasiparticles. This shows that although 
the bare particles are strongly interacting, the quasipar- 
ticles may not, and hints that the gas in the normal phase 
is probably behaving as a Fermi liquid, even very close 
to the resonance. 

In conclusion, we have built a theory which takes 
into account the many body exchange effects in the 
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quasiparticle-quasihole channel. This theory, contrary to 
the RPA. respects the Pauli principle and does not give 
spurious ground states. Inclusion of the exchange effects 
is therefore fundamental in obtaining the correct physics. 
We obtain a finite scattering amplitude as seen experi- 
mentally. We have also shown that a triplet superfluid 
is possible within the temperatures today achievable in 
cold atom traps and that triplet paring should be taken 
into account when discussing the properties of the sys- 
tem close to resonance. Furthermore, it seems possible 
in this formulation to derive the basic properties on the 
BEC side, although one should include properly the pres- 
ence of the bound state, which we have not. The strong 
agreement with experiments indicates that quasiparti- 
cles, not bare particles, are binding in the open channel. 
The good interpolation of the intermediate interacting 
region between the BCS and BEC sides is probably due 
to a careful account of the particle-hole contributions in 
the theory. Finally, we remark that this approach, since 
it is not restricted to the dilute gases, can be applied to 
other systems. 
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